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ABSTRACT 

We present high-resolution Chandra X-ray observations of PSR B0833 — 45, the 89 ms pulsar associ- 
ated with the Vela supernova remnant. We have acquired two observations separated by one month to 
search for changes in the pulsar and its environment following an extreme glitch in its rotation frequency. 
We find a well-resolved nebula with a toroidal morphology remarkably similar to that observed in the 
Crab Nebula, along with an axial Crab-like jet. Between the two observations, taken ~ 3 x lO^s and 
~ 3 X lO^s after the glitch, the flux from the pulsar is found to be steady to within 0.75%; the 3ct limit on 
the fractional increase in the pulsar's X-ray flux is ^ 10~^ of the inferred glitch energy. We use this limit 
to constrain parameters of glitch models and neutron star structure. We do find a significant increase 
in the flux of the nebula's outer arc; if associated with the glitch, the inferred propagation velocity is 
^ 0.7c, similar to that seen in the brightening of the Crab Nebula wisps. 

We propose an explanation for the X-ray structure of the Vela synchrotron nebula based on a model 
originally developed for the Crab Nebula. In this model, the bright X-ray arcs are the shocked termination 
of a relativistic equatorial pulsar wind that is contained within the surrounding kidney-bean shaped 
synchrotron nebula comprising the post-shock, but still relativistic, flow. In a departure from the Crab 
model, the magnetization parameter a of the Vela pulsar wind is allowed to be of order unity; this 
is consistent with the simplest MHD transport of magnetic field from the pulsar to the nebula, where 
B < A X 10"'' G. The inclination angle of the axis of the equatorial torus with respect to the line of 
sight is identical to that of the rotation axis of the pulsar as previously measured from the polarization 
of the radio pulse. The projection of the rotation axis on the sky may also be close to the direction of 
proper motion of the pulsar if previous radio measurements were confused by orthogonal-mode polarized 
components. We review effects that may enhance the probability of alignment between the spin axis and 
space velocity of a pulsar, and speculate that short-period, slowly moving pulsars are just the ones best- 
suited to producing synchrotron nebulae with such aligned structures. Previous interpretations of the 
compact Vela nebula as a bow-shock in a very weakly magnetized wind suffered from data of inadequate 
spatial resolution and less plausible physical assumptions. 

Subject headings: pulsars: general — pulsars: individual (PSR B0833 — 45) — X-rays: general — 
supernova remnant — stars: neutron 



1. INTRODUCTION 

Within two years of the discovery of radio pulses from 
CP1919-I-21, magnetized, rotating neutron stars were 
firmly established as the origin of these remarkable sig- 
nals. Furthermore, the steady increase in pulse period 
recorded for all sources provided an explanation for the 
pulsar power source: rotational kinetic energy. The de- 
tection of a decelerating 33 msec pulsar in the Crab Neb- 
ula solved the long-standing mystery of what powered this 
unique nebula: the spin-down rate of the Crab pulsar im- 
plied an energy loss rate, E 5x 10^^ erg s~^, more than 
enough to cover the radiation losses observed from radio 
to gamma ray frequencies. By 1974, the basic model of 
the electrodynamics of pulsar magnetospheres and their 
coupling to the surrounding synchrotron-emitting plasma 
was in place (Rees and Gunn 1974), although a detailed 
understanding of the processes involved continues to elude 
us (e.g., Arons 1998). 

After a decade of timing observations, it became clear 
that most pulsars were not defect-free clocks which sim- 
ply slowed smoothly as rotational energy was transformed 
into an electromagnetic outflow. Two distinct types of 
non-monotonic behavior were established: "timing noise" 



characterized by a stochastic wandering in pulse phase 
and/or frequency which appeared to afflict most pulsars 
(Helfand, Taylor, and Backus 1980; Cordes and Helfand 
1980; Cordes and Downs 1985; D'Amico et al. 1998), 
and "glitches" , an apparently instantaneous increase in the 
pulse frequency (a spin-up) accompanied by a simultane- 
ous change in the spin-down rate; these rare events were 
found to be most prevalent in young objects (Reichly and 
Downs 1971; Lyne 1996 and references therein). Thirty 
years after the first glitch in the Vela pulsar was recorded, 
a total of 65 events have been seen in 27 different pulsars 
(Lyne 1996; Wang et al. 2000). 

Glitches are a sudden fractional increases in the pulsar 
spin frequency with Sv/v w 10"^ to 6 x 10~^. No pulsar 
with a characteristic age of > 10^ yr has been observed to 
glitch more than once, but some young objects experience 
these events roughly annually. The best studied and most 
prolific in terms of large glitches is the first object in which 
a glitch was seen - the Vela pulsar. A dozen events have 
been recorded over the past three decades and daily mon- 
itoring continues. The largest event yet observed occurred 
in 2000 January (Jz^/z/ — 3.14 x 10~^) and provided the 
stimulus for the observations reported here. 
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The Vela Pulsar X-ray Nebula 



In this paper, wc report new results on PSR B0833 — 45 
based on obsc^rvations acquired with the Chandra High 
Resolution Camera. The data enable us for the first time 
to distinguish morphological details of the synchrotron 
nebula surrounding PSR B0833 — 45, and reveal a strik- 
ing picture of bilateral symmetry reminiscent of the loops 
and jets recently resolved in the Crab Nebula (Weisskopf 
et al. 2000). We offer an interpretation of the nebula's 
structure which requires an MHD wind with a high mag- 
netization parameter (unlike that inferred for the Crab). 
We also construct a high quality X-ray pulse profile and 
set tight upper limits on any change in the profile following 
the glitch, constraining models for the neutron star inte- 
rior. Finally, we demonstrate an apparent brightening in 
the Nebula a month after the spin- up event; whether this 
was stimulated by the glitch or is a phenomena akin to the 
flickering wisps in the Crab Nebula remains a question for 
future observations to answer. 

In section 2, we describe in some detail the analysis pro- 
cedures required to extract quantitative information from 
our HRC data, in part as a cautionary tale for other early 
users of this instrument. We then go on to explore the 
morphology of the pulsar's synchrotron nebula (§3), the 
soft X-ray pulse profile and limits on changes thereto (§4) , 
and a search for changes in the nebula following the glitch 
(§5). The Discussion (§6) begins by developing a model 
which accounts for the nebula's geometry with respect to 
the pulsar, as well as its anomalously low value of Lx/E; 
we then go on to derive constraints on glitch models from 
the temporal changes we see (and don't see) following the 
glitch. The final section (§7) summarizes our conclusions 
and assesses the prospects for future observations. 

2. OBSERVATIONS 

In response to an lAU Circular announcing a large Vela 
glitch on 2000 January 16.319 (Dodson et al. 2000), we 

submitted a Target of Opportunity request to the Chan- 
dra Observatory (Weisskopf, O'Dell, and van Speybroeck 
1996) to observe the pulsar as soon as practical, followed 
by a second observation roughly one month later in or- 
der to search for changes in the pulsar's flux, pulse profile, 
and/or surroimding nebula. The observations were carried 
out on 20 January 2000 and 21 February 2000, ~ 3.5 and 
~ 35 days after the glitch using the Chandra imaging High 
Resolution Camera (HRC-I; Murray et al. 1997). Integra- 
tion times of ~ 50 ksec were achieved in both observations. 

The HRC-I detector on-board Chandra is sensitive to X- 
rays over the 0.08 — 10.0 keV range, although essentially 
no energy information on the detected photons is avail- 
able. Photons are time-tagged with a nominal precision of 
15.6 /is; in this work, their arrival times were corrected to 
the solar system barycenter using a beta version of AXBARY. 
The data were collected during a portion of the orbit which 
avoided regions of high background contamination from 
the bright Earth and radiation belt passages; the second 
observation was, however, found to be partially contami- 
nated by particle activity, most likely of solar wind origin 
(see below). The pulsar was centered at the on-axis po- 
sition of the HRC where the point-spread function (PSF) 
has a minimum half-power diameter (the radius enclosing 
50% of total source counts) of ~ 0'.'5, which increases with 

^This problem does not affect early HRC-I observations obtained 



energy. Images were extracted centered on the pulsar and 
binned using the native HRC 0'.'13175 x 0'.'13175 pixel size 
into 1024 X 1024 pixel images (2.5' on a side). 

We began our analysis using event data calibrated by the 
initial processing and made available through the Chan- 
dra public archive. The first observation revealed several 
problems in the standard data sets and further problems 
were subsequently found during the analysis of the second 
observation. These problems affected both the spatial and 
timing analyses and had significant implications for the 
proper interpretation of the data. We alerted the HRC 
hardware and software teams to the instrument and data 
processing anomalies, and received considerable support 
in working through the problems. We document here the 
various artifacts discovered and the steps taken to correct 
for, or eliminate, them in our final data sets; our goals in 
doing so are 1) to allow for the replication of our results, 
and 2) to alert other early HRC users to problems they 
may encounter. In fact, we found it necessary to reprocess 
the data from Level 0.5 using custom scripts which incor- 
porated improved filtering and processing tools, making 
use of several beta versions of software provided by the 
HRC team. 

In the first observation we found evidence for a signif- 
icant "ghost image" which appeared as a spectacular jet- 
like feature emanating from the pulsar along the detector 
t;-axis. Examination by the instrument team found that 
the standard processing had failed to screen out all events 
flagged as instrumental. After filtering with a beta ver- 
sion of SCREEN_HRC with the mask parameter set to 32771, 
a much truncated jet-like feature was still apparent. In 
order to isolate any detector-centric artifacts, we obtained 
our second observation with a roll angle offset by 36° from 
the first. As discussed below, we were able to confirm the 
reality of the residual jet-like feature in the cleaned images. 

Independent of any filtering, the initial images showed 
the pulsar to be broader than the nominal PSF. To sep- 
arate the pulsar from any proximate nebular emission, 
we followed the same phase-resolved imaging analysis de- 
scribed in Gotthelf & Wang (2000) for their HRC obser- 
vation of PSR 0540-69, the 50 ms pulsar in the LMC. 
This separation, however, failed completely. By plotting 
the arrival times of the pidsar centroid, we observed that 
the sky coordinates of the pulsar wandered in a sinusoidal 
fashion with an amplitude of 0'.'3 and the periodicity of the 
programmed telescope dither. This accounted for the pul- 
sar's non-point-like appearance in the time-integrated im- 
age. Discussions with the Chandra attitude aspect team, 
however, showed a high-quality aspect reconstruction for 
the Vela observations. 

Further analysis by the HRC team revealed a systematic 
problem with one of the three anode preamplifiers which 
causes the coarse position algorithm to mis-place photon 
locations depending on the photon input position relative 
to the HRC tap gaps. This explains the apparent wan- 
dering of the pulsar centroid at the dither frequency: a 
fraction of the detected photons are displaced along the 
detector coordinates by a fixed amount. Indeed, we were 
able to ascertain that the apparent diffuse flux was pro- 
duced by faint echoes of the pulsar itself along the two 
orthogonal detector axes ^. 

prior to an increase in the instrument gain by a factor of two, such as 
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To eliminate the echoes, we initially used the bright pul- 
sar as a fiducial point to re-aspect the photons and thus 
take out the detector-induced wobble in a statistical sense. 
Subsequently, the instrument team made available a beta 
version of a code to identify and correct the mis-placed 
photons - HRC_EVTO_CORRECT ~ along with updated de- 
gap parameters (c/ul = 1.068; c/u2 = 0.0; cfvl = 1.045; 
cfv2 = 0.0) for use in HRC_PROCESS_EVENTS. This software, 
together with the new parameters, effectively eliminated 
the echo problem. The images produced by the two meth- 
ods are indistinguishable and the pulsar now matches the 
PSF to within its estimated uncertainty. 

In the second observation, we noted additional artifacts 
in the sky image resembling a rabbit-ear antenna extend- 
ing 20" from the pulsar in orthogonal directions along the 
detector axes with point-like sources at the ends. Tem- 
poral analysis of the region showed that the rabbit-ear 
counts occurred during a number of specific time inter- 
vals lasting tens of seconds. Examination of the mission 
timeline parameters showed that the occurrence of those 
events always followed the "AOFF_GAP" times by a few 
hundred seconds. Furthermore, data drop-outs were found 
for tens of seconds at the "AOFF_GAP" times and dur- 
ing the following intervals when the spurious counts were 
recorded. We wrote an algorithm to generate a now good- 
time-intervals file which eliminates these intervals based 
on the "AOFF.GAP" times. 

Two additional detector issues needed consideration 
when extracting accurate timing information: telemetry 
saturation and a hardware time-stamp mis- assignment. 
Although data obtained during the first observation dis- 
played nominal background levels, the second observation 
was plagued by intervals of telemetry saturation induced 
by high background levels; such occurrences can seriously 
affect timing studies by introducing spurious periods in 
the power spectrum aliased with the full buffer rate of 
~ 4 ms. We filtered out telemetry-saturated time inter- 
vals with the dead-time fraction criteria of DFT > 0.9. A 
further complication for precision timing was recently dis- 
covered by the HRC hardware team: the time-stamps for 
each event arc mis-assigned to the following event. Based 
on the VALID_EVT.COUNT count rate of ~ 500 cps, the 
average error in the assigned photon arrival time is 2 ms 
or a 2% phase error for the Vela pulsar; assuming roughly 
Poisson fluctuations in the HRC count rate over the obser- 
vation interval, the maximum error for any photon will be 
< 3 ms. Thus, with 25 phase bins across the 88 ms pulsar 
period, few, if any, photons have been misassigned and we 
have taken no mitigating action to correct this error. 

To compare directly the two observations, we repro- 
cessed both data sets starting from the Level 0.5 event files 
using identical methods and filter/screening/processing 
criteria, compensating for incorrect keyword values, pro- 
ducing correct GTI extensions, etc. This resulted in a total 
of 50.3 ks and 45.3 ks integration times for the first and 
second observation, respectively. Despite all the initial dis- 
crepancies and artifacts in the two observations, this repro- 
cessing produced effectively identical images, light curves 
and count rates. We are thus confident that we have elimi- 
nated all currently recognized instrumental artifacts in the 
final data sets upon which we base the analysis herein. 

those of the 50 ms pulsar PSR 0540-69 (see Gotthelf & Wang 2000). 



3. AN IMAGE OF THE VELA PULSAR 

A global view of the Vela pulsar and its environment as 
seen by the Chandra HRC is presented in Figure 1. The 
pulsar is embedded in a complex region of previously re- 
solved thermal X-ray emission from the Vela supernova 
remnant that is present throughout this image and ex- 
tends far beyond its boundaries. The X-ray jet noted by 
Markwardt and Ogelman (1995) is essentially overresolved 
in this image and extends far to the south of the image 
boundary; it is evident as a faint enhancement in the dif- 
fuse emission extending to the southeast and south of the 
bright pulsar nebula. 

The superb spatial resolution of the Chandra HRC pro- 
vides the first look at the structure of the synchrotron 
nebula in the immediate vicinity of PSR B0833 — 45; Fig- 
ure 2 shows an image constructed from the two observa- 
tions which have been centered on the pulsar and summed. 
The bright point source representing the pulsar has an ex- 
tent roughly consistent with the local PSF. Apparently 
emanating from the pulsar, towards the southeast, is a lin- 
ear, jet-like feature 10" in length. There is also evidence 
for a counter jet in the opposite direction. These jets have 
a position angle of 130° (measured East of North), and are 
aligned to within 8° ± 5° degrees with the pulsar's proper 
motion vector (Bailes et al. 1990; DeLuca, Mignani, and 
Caraveo 2000). 

Concentric with the pulsar is a diffuse outer arc of emis- 
sion perpendicular to the jet. This feature is roughly el- 
liptical in shape and subtends an angle of ~ 150 degrees 
as seen from the pulsar. Interior to this arc is an elliptical 
ring of emission with a curvature very similar to the outer 
arc. The pulsar, jet, and arcs are embedded in a extended 
nebula of faint diffuse emission which has been described 
as "kidney-bean" shaped (Markwardt and Ogelman 1998). 
The configuration of the jet feature relative to the nebula 
is reminiscent of the Chandra image of the Crab Nebula 
(Fig. 3; see Weisskopf et al. 2000). 

We determined the count rates by extracting counts 
from the various regions discussed above (see Figure 4). 
For each source region we carefully estimated the back- 
ground. For the diffuse emission, we determined the HRC 
detector background derived from an annulus 13.2" wide 
exterior to the kidney bean emission (r > 52.7"). We 
extracted counts from the pulsar using a circular aper- 
ture 2"64 in radius and estimated background from the 
surrounding annulus, 2"64 < r < 3" 43. Table 1 summa- 
rizes the properties of the individual components, includ- 
ing their estimated sizes and intensities. 

The best current measurements for the Vela pulsar and 
remnant place it at a distance of only 250 ± 30 pc (Cha, 
Sembach, and Banks 1999 and references therein); in all 
that follows we scale hy d = 250c?250 pc- 

Vela X, the ^ 100' diameter, flat-spectrum radio com- 
ponent near the center of the Vela remnant (Milne 1968; 
Bock et al. 1998) is generally regarded as the pulsar's 
radio synchrotron nebula. While soft X-rays from this re- 
gion are detected, they are primarily thermal in nature, 
and represent emission from the hot plasma which fills the 
entire remnant (Kahn et al. 1985; Lu and Aschenbach 
2000). Even the bright radio filament detected by Bieten- 
holz, Frail and Hankins (1991) shows no corresponding en- 
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Table 1 

HRC Spatial Components of PSR B0833 - 45 



PSR B0833-45 System 
Component" Shape and Size Counts s 



A) Pulsar (Obs 


1) 


point-like 


2.012 ±0.006 


(Obs 


2) 




1.997 ±0.007 


B) Nebula (Obs 


1) 


20" X 10" NE-SW 


2.720 ±0.007 


(Obs 


2) 




2.762 ±0.008 


C) Arc (Obs 


1) 




0.675 ± 0.004 


(Obs 


2) 




0.700 ± 0.004 


D) Jet (Obs 


1) 


10" long SE-NW 


0.037 ±0.001 


(Obs 


2) 




0.036 ±0.001 



"^See §3. 

''Count rates are background subtracted. 



hancement in our X-ray image, although more constrain- 
ing limits will be derivable from ACIS observations. 

The compact X-ray source near the pulsar was first rec- 
ognized by Kellogg et al. (1973) as having a harder spec- 
trum. Subsequent observations with increasing angular 
resolution (Harnden et al. 1985 and references therein; 
Ogelman, Finley, and Zimmermann 1993; Markwart and 
Ogelman 1998) localized the compact nebula to a region 
~ 2' in extent roughly centered on the pulsar. Ogelman et 
al. (1993) used the nominal PSF of the ROSAT PSPC 
and an a,d hoc model for the surface brightness of the 
diffuse emission to attempt a deconvolution of the pul- 
sar and its nebula and to obtain spectral fits to the two 
components. They found that a blackbody efi^ective tem- 
perature of 0.15 keV adequately characterized the point 
source, while the extended emission exhibited a power law 
spectrum with a photon index of ^ 2.0; a column den- 
sity of Nh = 1 X 10^° cm~^ is marginally consistent with 
both components. Markwardt and Ogelman (1998) subse- 
quently revised the division of the flux between the point 
source and nebula based on ROSAT HRI observations. Se- 
ward et al. (2000) attempted to isolate the pulsar emission 
temporally and found somewhat lower blackbody temper- 
atures. The RXTE observations of Gurkan et al. (2000) 
found a similar power law index for the nebula, but a much 
higher normalization; while this could indicate diffuse syn- 
chrotron X-ray emission from a larger area (given their 
one-degree field of view), it could also result from back- 
ground modeling problems, since Vela is a weak source for 
RXTE. 

As described above, om- HRC image allows us to sep- 
arate cleanly the pulsar from the nebular emission. We 
find count rates for the pulsar and the nebula minus 
the pulsar (within a radius of 50") completely consistent 
with those of Markwardt and Ogelman (1998) using the 
Ogelman et al. (1993) spectral parameters {kT = 0.15 
keV; Nh = 1 x 10^° cm^^), confirming these as a uscfid 
characterization of the observed fiux. This leads to an un- 
absorbed, bolometric luminosity for the pulsar blackbody 



emission of w 1.5 x lO'^^d^r^Q erg s^^. Note, however, that 
this luminosity is inconsistent with the adopted temper- 
ature (T = 1.7 x lO^K) and distance (250 pc) for a uni- 
formly radiating blackbody with a radius of 10 km; even 
for the minimum neutron star radius consistent with rea- 
sonable equations of state (i? ~ 7 km), the derived Lx is 
too high by a factor of 20. Lowering T to 8.5 x lO^K as 
advocated by Seward et al. (2000) yields self- consistent 
values for Lx , R, T, and d, and raises the intrinsic lumi- 
nosity by ^ 30%. Adjusting the distance, and including 
such effects as non-grey opacity in the neutron star atmo- 
sphere and a non-uniform temperature distribution over 
the surface will also aS^ect the calculated luminosity. For 
the purpose of a comparison with glitch models (§6), we 
adopt T = 1.0 X 10*^' K. 

The integrated luminosity of the whole nebula (r < 
52.7" minus the pulsar) in the 0.1 — 10 keV band is 
3.5 X lO'^^dgso erg s~^, corresponding to 4.9 x 10~^ of the 
pulsar's spin down luminosity. This ratio of Lneb /E is sig- 
nificantly lower than that for any other pulsar, and is a 
major constraint on models for coupling the pulsar wind 
to the nebula (see §6). 

4. THE X-RAY PULSE PROFILE 

After many unsuccessful searches, Ogelman et al. (1993) 
were the first to detect X-ray pulsations from the Vela pul- 
sar using the ROSAT PSPC. Their observations revealed 
a complex profile, not obviously related to the pulse pro- 
files previously recorded at radio, optical, and gamma-ray 
wavelengths. The observed pulsed fraction of 4.4 ± 1.1% 
was diluted by the inability of the PSPC to resolve the pul- 
sar from the surrounding nebula; using the approximate 
model described above, the authors estimated a soft X-ray 
pulsed fraction of 11%. Seward et al. (2000) constructed 
a higher signal-to-noise profile by combining six ROSAT 
HRC observations; they estimated a pulsed fraction of 12% 
divided between a broad component (8%) and two narrow 
peaks (4%). Strickman et al. (1999) and Gurkan et al. 
(2000) have recently derived 2-30 keV profiles based on 
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RXTE observations of Vela. Strickman et al. illustrate 
a trend in which the component separation of the main 
pulse increases with energy. 

We have determined the X-ray pulse period for our two 
observations and compared them to the radio ephcmcris 
(Don Backer, personal communication). We began by con- 
structing a pcriodigram around a narrow range of periods 
centered on the expected period ±0.1 ms, sampled in in- 
crements of 0.05 X P^/T, where T is the observation du- 
ration, and P is the test period. For each trial period, 
we folded photons extracted from a r = 2'.' 64 aperture 
centered on the pulsar position using 25 phase bins and 
computed the of the resultant profile. We find a highly 
significant signal (> Scr) at P = 89.32842(5) ms at epoch 
51563.314043 M.JD (TBD) and P = 89.32876(6) ms at 
epoch 51595.370251 MJD (TBD), completely consistent 
with the observed radio periods. The uncertainty was esti- 
mated according to the method of Leahy (1987). For each 
fold, we adoped the period derivative determined from the 
radio ephcmcris. Our measurements of the period are con- 
sistent (within the errors) with the radio prediction. 

In order to compare the pulse profiles for the two ob- 
servations we used the phase-connected radio ephemeris 
to fold and align them. We present the sum and differ- 
ence profiles in Figure 5. This phase alignment is com- 
pletely consistent with one we computed empirically by 
cross-correlating the two profiles; this suggests that the 
Chandra clock is stable to a few ms over a month. Wc 
can also derive the absolute radio to X-ray phase offset 
if we assume that the absolute Chandra time assignment 
is accurate (the calibration of this quantity has not yet 
been finalized). The phase of the radio peak relative to 
the X-ray profile is indicated in Figure 5. 

The ~ 200,000 counts, uncontaminated by nebular 
emission, provides us with the highest signal-to-noise X- 
ray pulse profile for Vela yet reported (Figure 5). Greater 
than 99% of the counts from the blackbody component de- 
tected by the HRC fall in the 0.1 - 2.4 kcV ROSAT band. 
We compute a pulsed fraction by integrating the counts in 
the light curve above the lowest point and dividing by the 
total counts within a radius of r = 2.64" and subtracting 
the small amount of nebular background in this extraction 
radius (see Table 1). Our value for the pulsed fraction is 
7.1 ± 1.1%; the quoted error is dominated by the Pois- 
son uncertainty in the number of counts recorded in the 
bin representing the light curve minimum. This value is 
somewhat lower than those cited above, but has the advan- 
tage of utilizing a direct measure of the total point-source 
contribution with subarcsecond resolution. We measure a 
separation between the two peaks of the main component 
of 54> ~ 0.325, consistent with a linear extrapolation of the 
energy-dependence of this quantity reported by Strickman 
et al. (1999). 

The background-corrected pulsar count rates were found 
to be 2.012 ± 0.006 and 1.997 ± 0.007 c s-\ respectively, 
for the first and second epochs; the overall comit-rate is 
constant to within 0.75% (the second observation is 1. So- 
fainter than the first). Thus, the 3(T limit on any increase 
in the pulsar luminosity in response to energy input from 
the glitch is < 1.2 x lO^o erg s'^ or AT ~ o'.2%, 35 days 
(3 X 10^ s) after the event^. The lower half of Figure 5 

^This estimate includes a ~ 20% correction for the fraction of the 



shows the difference between the two observations as a 
func;tion of pulse phase, where the second dataset has sim- 
ply been scaled by the ratio of the total integration times; 
no single bin has a discrepancy exceeding 1.5 sigma. This 
constancy in both the pulsed luminosity and pulse profile 
set interesting constraints on the glitch mechanism (see 
§6.4). 

5. CHANGES IN THE NEBULA 

While the primary energy release from a glitch must be 
within (or on the surface of) the neutron star, the response 
of the star's magnetosphere could result in the release of 
energy to the synchrotron nebula, triggering changes in its 
morphology and/or brightness. Even without the stimulus 
of a glitch, the optical wisps of the Crab Nebula near the 
pulsar have been shown to change on timescales of weeks, 
presumably in response to instabilities in the relativistic 
wind from the pulsar (Hester et al. 1995). Greiveldinger 
and Aschenbach (1999) have also reported changes in the 
X-ray surface brightness of the Crab Nebula on larger 
scales and somewhat longer timescales. Thus, we have 
examined our two images of the Vela nebula carefully in a 
search for surface brightness fluctuations. 

We examined the count rate in the kidney-bean re- 
gion bracketed by the outer-arc and an inner circle with 
r = 1'.'32 centered on the pulsar. No significant change 
was observed between the two observations. Similarly, no 
measurable change in the count rate associated with the 
jet-like feature was found. Comparisons of other regions 
defined in Table 1 also showed no change, with the notable 
exception of the outer arc itself which appeared to increase 
in brightness by ^ 5% in the second observation. 

To investigate this further, we examined regions con- 
gruent with the morphology of the nebula by constructing 
radial bins which are elliptical in shape; the ratio of the 
semi-major to semi-minor axes of the ellipse is 1.76, and 
the elliptical annuli are oriented at a position angle of 50° 
(east of north). As Figure 6 shows, the sector of the radial 
profile encompassing the bright northwestern arc exhibits 
a 7.8cr excess between semi-minor axis radii of 13.5" and 
18.0" in the sense that the source is brighter in the second 
observation. No other sectors or radii show any significant 
changes. In Figure 7, we display the azimuthal profile of 
the whole nebula in the elliptical ring 4.5" wide centered 
on these radii. The residuals are positive (the second ob- 
servation is brighter) throughout the range 250° to 30°; 
the excess is significant at the 7.6a level and represents a 
brightening of 5.3%. The excess energy being radiated in 
the HRC band amounts to ^ 3 x 10^° erg s~-^. For the 
assumed geometry discussed in §6.1, the outer arc lies at a 
distance of 1.05 x 10^^ cm from the pulsar, requiring signal 
propagation at ^ 0.7c if the impetus for the brightening 
originated from the pulsar at the time of the glitch. 

6. DISCUSSION 

6.1. Geometry and Kinematics of the Nebular Structure 

Figure 8 shows our proposed model for the Vela X-ray 
nebula. We assume that the two prominent arc-like fea- 
tures lie along circular rings highlighting shocks in which 
the energy of an outflowing equatorial wind is dissipated 
to become the soiirce of synchrotron emission for the com- 

metric luminosity lying below the HRC band. 
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pact nebula extending to the boundary of the "bean" . One 
reason that the arcs are not complete rings might be that 
the emission is from outflowing particles which Doppler 
boost their emission in the forward direction. 

This is essentially the picture for the similar arcs sur- 
rounding the Crab pulsar first suggested by Aschenbach 
and Brinkmann (1975) and later elaborated by other au- 
thors (Arons et al. 1998 and references therein). The 
main difference is that the dark cavity which contains the 
imshocked pulsar wind in the Kennel & Coroniti (1984a) 
model of the Crab is small compared with the volume of 
the Crab Nebula, while the radius of the Vela nebula (the 
bean), is barely twice as large as its pulsar wind cavity 
(see Figure 4). We assume that the two rings straddle 
the equator symmetrically, and suppose that the deficit 
of emission exactly in the equatorial plane is related to 
the fact that this is where the direction of a toroidally 
wrapped magnetic field changes sign; i.e., the field may 
vanish there. The semimajor axis of the ring a = 25.7", 
and the ratio a/b = 1.67 specifies the angle that the axis 
of the torus (i.e., the rotation axis of the pulsar) makes 
with the line of sight, ^ = cos~^(&/a) = 53°. 2. The angle 
^'o = 130° is the position angle of the axis of the torus on 
the plane of the sky, defined according to convention as 
the angle measin-ed to the east from north. The direction 
of rotation (the sign of f2) is arbitrary. 

The projected separation of the two rings is measured 
as s = 17.7". The half opening angle of the wind 9 is 
then given by tan6' = s/(2asinC) = 0.43(61 = 23.°3), 
and the radius of the shock is rg — ad /cos 9, where d 
is the distance to the pulsar. For d = 250 pc we find 
rs = 1.05 X 10^'^ cm. 

The rotating vector model of pulsar polarization (Rad- 
hakrishnan & Cooke 1969) is commonly used to derive 
information about the geometry of the pulsar magnetic 
inclination and viewing angles. The angles ( and can 
in principle be evaluated independently using information 
derived from polarization measurements of the radio pulse. 
In particular, the swing in position angle \l/(i) of linear po- 
larization across the pulse is very sensitive to (, the angle 
between the line of sight and the rotation axis. The an- 
gle a between the magnetic axis and the rotation axis is 
much more difl&cult to measure unless a » i.e., unless 
the line of sight passes near the center of the polar cap. 
Accordingly, a is often assumed while ( is fitted. For ex- 
ample, Krishnamohan & Downs (1983) assume a = 60° 
in their model for Vela; this is consistent with the value 
a = 65° derived by Romani and Yadigaroglu (1995) from a 
fit of their geometric gamma-ray emission models to Vela's 
pulse profile. When an interpulse is observed, a is often 
inferred to be 90° (that is, both polar caps are visible in 
this case). With these definitions (see Figure (8) or Figure 
(13) of Krishnamohan & Downs (1983) for the geometry), 



tan(*(i) - *o) 



sin (j){t) 



cot a sin — cos^ cos^(t) 



(1) 



Here (f>(t) is the longitude of the emitting region, which 

increases linearly with time, and 5*0 is the position angle 
of the rotation axis of the pulsar projected on the sky as 
in Figure (8). 

In the context of the rotating vector model, ^'o is iden- 
tical to the position angle of polarization ^ at the peak of 



the pulse where the magnetic dipole axis crosses the rota- 
tion axis {(t){t) = in Equation (1)). Because the emission 
mechanism is thought to be curvature radiation from par- 
ticles moving along magnetic field lines, the electric vector 
is tangent to those field lines, rather than perpendicular 
to them as is the case with synchrotron radiation. While 
this measurement is in principle straightforward, in prac- 
tice it is not routinely accomplished. Observations at two 
or more frequencies are needed to determine (and to cor- 
rect for) the interstellar rotation measure, and to demon- 
strate that the intrinsic polarization is in fact frequency- 
independent. Another complication is that a pulse is of- 
ten composed of several identifiable components, some of 
which can be polarized in the orthogonal mode (a result of 
propagation effects in the magnetosphere) obscuring the 
"true" polarization. Furthermore, the pulse itself might 
not even contain an identifiable core component, being 
composed instead of emission from random patches within 
a cone (e.g., Lyne and Manchester 1988; Manchester 1995; 
Deshpande and Rankin 1999). Accordingly, measurements 
of \E'o are rarely attempted, measurements of a are rarely 
trusted, and measurements of (J are rarely questioned. 

In fact, there are several determinations of for the 
Vela pulsar that are not in particularly good agreement 
with each other; we review a representative subset here. 
The original value of Radhakrishnan & Cooke (1969) is 
^0 = 47° with an uncertainty of w 5°. Hamilton et al. 
(1977) made measurements over several years, all of which 
are consistent with 5'o = 64° ± l.°5. A detailed decompo- 
sition into four separate pulse components was performed 
by Krishnamohan & Downs (1983), in which they con- 
cluded that one of the components was polarized in the 
mode orthogonal to the other three. However, they did 
not attempt an absolute measurement of the angle \E'o. Bi- 
etenholz et al. (1991) measured ^0 = 35° using the VLA 
(although this measurement may not be directly compa- 
rable, in that it represents a mean value weighted by the 
degree of linear polarization rather than the value at the 
center of symmetry of the position angle curve). Thus, 
while the published values of ^"0 differ by as much as 30°, 
it appears that none is even close to being aligned with the 
axis of the nebula, and that all are roughly perpendicular 
to it. 

Interestingly, the model of Krishnamohan & Downs 
(1983) produces a precise (albeit model-dependent) value 
of the angle between the rotation axis and the line of sight, 
C = 55.°57± 0.°15. This value agrees well with the inclina- 
tion angle of our postulated equatorial wind torus to the 
line of sight derived by fitting an ellipse to the shape of 
the X-ray features: ( = cos~^(6/a) = 53°. 2. This striking 
coincidence gives us courage to pursue the basic physics of 
the equatorial wind model using the geometry of Figure 8, 
and even to be so bold as to suggest that all of the radio 
determinations of for Vela are incorrect by 90° because 
of incorrect mode identification {i.e., perhaps three out of 
four of the pulse components are actually polarized in the 
orthogonal mode). In this case, = 130° (as inferred 
from the orientation of the X-ray torus), and can be iden- 
tified with the projected direction of the pulsar rotation 
axis. Speculations about the true orientation of \l/o in pul- 
sars go back to Tademaru (1977), who first discussed the 
possible alignment between spin axis and proper motion 
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in the context of the radiation rocket hypothesis (Harrison 
& Tademaru 1975). 

6.2. Implications of the Proper Motion 

For both the Crab and Vela pulsars, the direction of 
proper motion (transverse velocity Vt) is strikingly close to 
the projected X-ray symmetry axis of the inferred equa- 
torial wind and polar jet-like structures. The proper mo- 
tion of the Vela pulsar has been measured with compara- 
ble accuracy using both radio intcrferomctry (Bailes ct al. 
1990) and images from HST (De Luca et al. 2000). The 
resulting mean value 0."056 dz 0."004 yr~^ at a position 
angle of 302° ± 4° is within 8° of the axis of the neb- 
ula (^-o + 180° = 310°). The transverse velocity vi = 65 
km s^^ at = 250 pc. The proper motion of the Crab pul- 
sar is 0."018 ± 0."003 yr-i at a position angle of 292° ± 10° 
(Caraveo & Mignani 1999), which corresponds to Vt = 123 
km s~^ at d = 2000 pc. The axis of the Crab's toroidal 
optical and X-ray structure is 299° (Hester et al. 1995), 
only 7° from the direction of proper motion. The proba- 
bility that two such close alignments will occur by chance 
when drawn from a pair of uncorrelated distributions is 
0.7%. We also note that both Vela and the Crab are rather 
slow moving compared to typical young pulsars; Lyne & 
Lorimer (1994) found a mean velocity for young pulsars of 
between 400 and 500 km s~^. 

If these relationships are not a coincidence, then they 
may be understandable in terms of the scenario proposed 
by Spruit & Phinney (1998) who suggested that the rota- 
tion axes and space velocities of pulsars could be connected 
through the nature of the "kicks" given to neutron stars at 
birth. Spruit & Phinney argue that the rotation rate of the 
progenitor stellar core is too slow in the few years before 
the formation of the neutron star for pulsar spin periods 
to be explained by simple conservation of angular momen- 
tum during core collapse. Instead, it is likely that the same 
asymmetric kicks (whatever the cause) that are responsible 
for the space velocities of pulsars, are also the dominant 
contributors to their initial spin rates. If neutron stars 
acquire their velocities from a single momentum impulse, 
then their rotation axes should be perpendicular to their 
space velocities. If, however, they receive many random, 
independently located impulses over time, as might result 
from convection which leads to anisotropic neutrino trans- 
port or anisotropic fallback, then their velocities and spins 
should be uncorrelated in direction. However, if those mul- 
tiple thrusts are not short in duration relative to the re- 
sulting rotation period, it is possible that kicks applied 
perpendicular to the rotation axis will average out, while 
those that are along the rotation axis will accumulate. In 
the latter case, particularly germane for short rotation pe- 
riods, the space velocity will be preferentially aligned with 
the rotation axis. This is actually the scenario preferred by 
Spruit & Phinney, for which they appeal to long duration 
(several-second) thrusts that could result from the effect of 
parity violation in neutrino scattering in a magnetic field. 

Thus, in the context of the above models, we speculate 
that the Crab and Vela pulsars have relatively low space 
velocities because the components of their kicks perpen- 
dicular to their rotation axes were averaged out and, as a 
result, their final space velocities were aligned closely with 
their spin axes (cf. Lai, Chernoff, and Cordes 2001). We 
also note that alignment of the spin axis and proper motion 



is a natural consequence of the Harrison and Tademaru 
(1975) photon rocket acceleration mechanism. With the 
recent revision of Lai et al. (2001), an initial spin period 
as long as 6 ms suffices to account for the measured trans- 
verse component of the velocity in the maximum acceler- 
ation case. Since pulsars that rotate most rapidly at birth 
are also the ones most capable of powciring synchrotron 
nebulae, either scenario might argue for a stronger than 
average correlation of the axes of such nebulae with the 
proper motion directions of their parent pulsars. 

Before the toroidal arcs in the Vela nebula were resolved 
by Chandra, Markwardt & Ogleman (1998) interpreted the 
overall shape of the nebula as seen by the ROSAT HRI as 
being determined by the space velocity of the pulsar. In 
particular, they noted that the outline of the nebula, which 
is dominated by the bean shape, resembles a bow shock 
whose symmetry axis at position angle 295° is identical to 
the direction of the radio proper motion (297°). However, 
the rather uniform and gently curved outline of this struc- 
ture could only be reconciled with the sharper, asymmetric 
curve expected of a bow shock if the space velocity were 
nearly along the line of sight. This requirement, coupled 
with the large absolute velocity needed for the pulsar to 
exceed the speed of sound in the surrounding supernova 
remnant forced Markwardt & Ogleman to conclude that 
the velocity vector is less than 22° from the line of sight. 
This notion of the compact Vela nebula as a bow shock also 
led Chevalier (2000) to a considerably different model of 
its physics. As we shall argue below, the Chandra observa- 
tions do not support such a bow-shock interpretation, but 
instead favor a physical model in which the entire struc- 
ture is a synchrotron nebula similar in physics to the Crab, 
but with an interesting difference in one of its parameters. 

6.3. A Physical Model of the Nebula 

In the basic Kennel & Coroniti (1984a, b) model of the 
Crab Nebula, a relativistic pulsar wind terminates in an 
MHD shock, which produces the nonthermal distribution 
of particles and post-shock magnetic field that comprise 
the synchrotron nebula. Although the pulsar wind is as- 
sumed to carry the entire spin-down luminosity of the pul- 
sar, specific wind parameters such as the particle velocity 
and the fraction of the power carried in magnetic fields 
are not known a priori. Rather, they are inferred by us- 
ing the results of the shock jump conditions to model the 
spectrum of the Crab Nebula, and also to match the ob- 
served radii of the MHD shock and the outer boundary 
of the Nebula. It is necessary to adopt outer boundary 
conditions; a natural one is to require the final velocity of 
the flow to match the observed expansion velocity of the 
Nebula, although it is not clear how this outer boundary 
condition is communicated to the inner MHD shock which 
is a factor of 20 smaller in radius. 

A peculiar result of the Kennel & Coroniti model is 
that the wind magnetization parameter a is required to 
be w 0.003 in the Crab. That is, the fraction of power 
carried in B field is much less than 1%. Such a small frac- 
tion is required in order that sufficient compression occurs 
in the shock to convert the bulk flow energy into random 
energy of the particles so that they can radiate the ob- 
served synchrotron luminosity. Highly magnetized shocks 
produce less radiation because there is little energy dis- 
sipation and, for the Crab, would supply an insufficient 
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number of X-ray emitting particles. Furthermore, highly 
magnetized shocks are weak because all of the energy dis- 
sipation allowed by the jump conditions is used in making 
the small increase in B field needed to conserve magnetic 
flux. The post-shock flow velocity is still relativistic. 

The reason that such a small magnetization is difficult 
to understand is that the pulsar magnetic field energy car- 
ried out to the radius of the shock in the simplest MHD 
wind should be of the same order of magnitude as the spin- 
down power, as the following argument shows. The total 
wind energy flux at the shock is 

inn _ fBiRHi*\ j_ _ Bl n^ 

while the transported pulsar magnetic field Bg at the lo- 
cation of the shock is 

n (rV fnA _ B.R^n^ 

where ric is the radius of the light cylinder defined as 
Tic = c/n and R is the neutron star radius at which the 
magnetic field strength is Bp. Therefore, the magnetic 
energy flux at is 

^2 Bl i?6 1]4 

Stt ~ 327reV2 ' ^ ' 

comparable to the value in Equation (2). While various so- 
lutions to this paradox for the Crab have been proposed, 
we argue here that the dimensions and spectrum of the 
Vela synchrotron nebula are in much better accord with 
£7 ~ 1. 

A basic application of the pulsar wind model to the 
Vela synchrotron nebula was made by de Jager, Harding, 
& Strickman (1996) in conjunction with their detection 
of Vela with the Oriented Scintillation Spectrometer Ex- 
periment (OSSE) on the Compton Gamma-Ray Observa- 
tory. We summarize their conclusions here, de Jager et 
al. noted that the unpulsed part of the hard X-ray spec- 
trum extends with power-law photon index T = 1.73 up 
to 0.4 MeV, which implies that the synchrotron nebula ra- 
diates - 2 X 10''''(£;max/0.4MeV)°-2'^ ergs s-\ or only ~ 
3 X 10""^ of the pulsar spin-down power. The absence of an 
observed spectral break limits the residence time Tr of the 
electrons in the nebula that radiate in this energy range to 
less than their synchrotron lifetime t,, = 5.1 x 10*/ {j B'^) s. 
Since the electrons radiating at the highest observed en- 
ergy -Bmax have 7 = [ZTrmcEu-iax/iheB)]^/^, the upper 

2/3 

limit on the nebular magnetic field is _B < 20 r,- G. 
To make use of this limit, de Jager et al. assumed that 
the residence time would be r,. w r/cs where the nebula 
outer radius r ~ 2 x 10^^ cm, and Cg = c/VS, the ve- 
locity of sound in a relativistic plasma. For these values, 
B < A X 10~'^ G. This requirement is compatible with a 
strongly magnetized pulsar wind for which little change 
in B occurs across the shock, while the post-shock field 
continues to decline inversely proportional to the distance 
from the pulsar. Indeed, equation (3) predicts a pre-shock 
field of ~ 1.5 X 10~^ G, consistent with the limit derived 
above; thus, the wind may remain relativistic across the 



nebula, which extends a factor of 2 in radius beyond the 
shock. 

Such a synchrotron nebula is in approximate pressure 
balance with the surrounding supernova remnant. Mark- 
wart & Ogelman (1997) foimd by fitting a two-temperature 
thermal model to the ASCA spectrum of the inner rem- 
nant that the thermal pressure is « 8.5 x 10~^° erg cm~^. 
This compares well with the pulsar wind pressure at r, 
E / {4:TT r"^ c) = 5.2 x 10"^^° erg cm^'\ The entire compact 
X-ray nebula, then, is consistent with being powered by 
a strongly magnetized pulsar wind shock whose still rel- 
ativistic downstream flow is conflned by the Vela SNR. 
The modest production of synchrotron electrons in such 
a shock naturally explains the very low value of Lx / E 
observed. 

In summary, the match between the nebular field up- 
per limit derived from the OSSE data and the value found 
from equation 3 using our measured value of r^, and the 
fact that the pressure corresponding to this field strength 
matches the confining thermal pressure of the X-ray gas 
leads us to conclude that our high-magnetization model 
of the Vela nebula is both self-consistent and plausible. 
Thus, in a sense. Vela may be a more natural realization 
of the Kennel & Coriniti model than is the Crab, for which 
the model was created. 

While this appears to be a quite satisfactory model, cer- 
tain details are subject to additional constraints. First, 
since the radiation from the nebula is so inefficient, en- 
ergetic electrons must be able to escape to much larger 
distance scales before losing all of their energy to syn- 
chrotron radiation. Indeed, the radio luminosity of the 
100-arcminute Vela X region, ~ 8 x lO''^ erg s~^, could 
be one manifestation of the escaping electrons. Second, 
there is a natural upper energy to the synchrotron spec- 
trum when the electron gyroradius rg exceeds the radius 
of the nebula. Since 

(D \ -3/2 / pj \ -^/■^ 
TF^j (looM^j '^^'^ 

this is not a restrictive limit. Thus, we find the de Jager 

et al (1996) description of the Vela synchrotron nebula 
basically in accord with the Chandra observations. 

An alternative picture was proposed by Chevalier (2000) 
based on the bow-shock interpretation of Markwardt & 
Ogleman (1998), and a simplified version of the Kennel 
& Coroniti model. If the Vela synchrotron nebula is en- 
ergized by a shock between a relativistic pulsar wind and 
the surrounding supernova remnant, then the bow shock 
travels at the velocity of the pulsar Vp, which necessarily 
exceeds the sound speed in the hot confining medium. In 
the first place, such a large pulsar velocity is hardly likely, 
since the thermal sound speed y/dP/dp « 875 km s~^ 
in the Vela SNR according to the ASCA spectral analysis 
of Markwardt & Ogleman (1997). Markwardt & Ogleman 
(1998) assumed that Vp > 260 km s~^ is sufficient. 

An additional consequence of this scenario, however, is 
that the residence time of the emitting particles in the 
nebula is much longer than in the de Jager et al. model, 
Tr ^ r/vp instead of ~ \/3r/c. Consequently, Cheva- 
lier (2000) was forced to assume Tr ^ 10^ yr, requiring 
an extremely small magnetization parameter, a < 10~*, 
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in order that the model's radiated luminosity not exceed 
the observed X-ray luminosity. We consider that the phys- 
ical difficulties of the bow-shock interpretation, in conjunc- 
tion with the new Chandra evidence that the X-ray mor- 
phology is dominated by a pair of arcs resembling similar 
toroidal structures in the Crab Nebula, strongly disfavor 
such a model. Chevalier's revised estimate of a is 0.06 
after adopting our assumption of relativistic post-shock 
flow. Determining whether a is actually of order unity or 
not will require more detailed modeling 

While our model provides a plausible explanation for 
the torroidal arcs and outer nebula, it says nothing about 
the other striking feature of the image, the jet and coun- 
terjet. The jet is 10" long, giving it a deprojected length 
of 4.1 X lO""^^ cm. At w = 65 km s~^, it has taken the 
pulsar ~ 200 yr to travel this distance. The synchrotron 
lifetime of electrons producing 1 keV emission ranges from 
5 yr to 40 yr for the fields of (1 — 4) x 10"" G discussed 
above. Thus, the jet is not simply a wake, but must be 
supplied with particles from the pulsar (as the existence of 
the counterjet also suggests). The luminosity required is 
modest: Lx ~ 5 x 10^° erg s~^, roughly 1% of the nebular 
luminosity and < 10~^E. While we have no scenario to 
propose, the existence of such features in both Vela and 
the Crab suggests understanding their origin may prove 
useful in modeling particle flow from young pulsars. 

6.4. Thermal Emission Constraints on the Neutron Star 

Interior 

Several models have been advanced to explain the sud- 
den apparent change in the moment of inertia of a glitching 
neutron star. Originally, starquakes, resulting from the re- 
lease of strain in the stellar crust induced by the change 
in the equilibrium ellipticity of the star as it slows, were 
invoked (Ruderman 1969). But the magnitude and fre- 
quency of the Vela glitches could not be explained by this 
model, and a picture involving the sudden unpinning from 
the inner crust of superfluid vortices in the core of the 
star became the dominant paradigm (Anderson and Itoh 
1975). Observations of the relaxation of the star back to- 
ward its original spin-down rate suggest that ^ 1% of the 
star's mass is involved in the event (Alpar et al. 1988; 
Ruderman, Zhu, and Chen 1998), implying a total energy 
release of ~ 10"^^ ergs. 

The fate of this energy is unclear and predictions con- 
cerning the observable consequences vary widely. The 
timescale for energy deposited at the base of the crust to 
diffuse outward, the fraction of the surface area whose tem- 
perature will be affected, and the secondary effects, such 
as the rearrangement of the surface magnetic fleld which 
could dump energy into the surrounding synchrotron neb- 
ula, are all uncertain by one or more orders of magnitude. 

Our stringent upper limit on a change in the X-ray flux 
from the neutron star within the 35 days following the 
glitch allows us to begin setting meaningful constraints on 
the parameters of the neutron star and the glitch. Se- 
ward et al. (2000) provide a concise introduction to the 
published models for the thermal response of the stellar 
surface to a glitch generated by the sudden unpinning of 
superfluid vortex lines deep in the star (Van Riper, Ep- 
stein, and Miller 1991; Chong and Cheng 1994; Hirano et 
al. 1997, and Cheng, Li, and Suen 1998), and we need 
not repeat it here. We follow their approach in deriving 



parameter limits from the models. 

The allowable change in the pulsar's flux derived in §4 
corresponds to a fractional change in the surface temper- 
ature of < 0.2%. For timescales of ~ 30 days, we are 
primarily sensitive to stars with small radii (i? < 14 km) 
which correspond to soft or moderate equations of state. 
Using Figure 2 of Van Riper et al. (1991), we can set a 
limit of E glitch < 10^^ ergs independent of the depth of 
occurrence within the inner crust. For depths correspond- 
ing to local densities p < 10^"^ gm cm~^, Egutch < 3 x 10'*"'^ 
ergs, and for shallow events (p ~ 10^^ gm cm~^), the glitch 
energy must be less than 10"^^ ergs. For the softest equa- 
tion of state used by Hirano et al. (1997) corresponding to 
a 1.4M0 star with a radius of 11 km, glitch depths shal- 
lower than 10"'^^ gm cm^^ require an energy deposition of 
less than 10''^ ergs. With observations of similar sensitiv- 
ity ^ 300 and '~ 3000 days after the event, we could rule 
out Egiitch ~ 10^^ erg for all equations-of-state and glitch- 
depth combinations, and require Eguich < 10^^ ergs for 
soft and moderate equations of state for depths p < lO^'^"'' 
gm cm~^. Note that on the longest timescales (appropri- 
ate for deep glitches in stars with very stiff equations of 
state), Vela may be an inappropriate target for constrain- 
ing glitch parameters, since another glitch may well have 
occurred before the thermal pulse has peaked; indeed, if 
even 10% of the glitch energy appears as surface thermal 
emission, the total X-ray luminosity can be powered by 
events with < EgUtch >~ 2 x 10*^ ergs. 

7. SUMMARY AND CONCLUSIONS 

We have presented a high-resolution X-ray image of the 
Vela pulsar revealing a highly structured surrounding neb- 
ula. We interpret the nebula's morphology in the context 
of the shocked MHD wind model developed by Kennel and 
Coroniti (1984a, b) for the Crab Nebula, and find that the 
Vela nebula allows a large magnetization parameter, pos- 
sibly of order unity. This picture also provides a natural 
explanation for the low Lx/E of the Vela nebula. We 
speculate that the alignments of the symmetry axes of the 
Crab and Vela nebiflae with the proper motion vectors 
of their respective pulsars should be expected preferen- 
tially in rapidly spinning young pulsars with surrounding 
X-ray synchrotron nebulae if the causal connection be- 
tween spin and proper motion suggested by Spruit and 
Phinney (1998) is correct. 

Our two observations, centered 3.5 and 35 days after 
the largest glitch yet recorded from the pulsar, allow us 
to set significant limits on changes in the pulse profile and 
stellar luminosity which can be used to constrain glitch 
model parameters. We find that, for soft and moderate 
equations of state, the glitch energy must be < 10"*^ ergs; 
an additional observation a year following the event will 
substantially tighten this constraint. An apparent change 
in the nebula surface brightness between the two obser- 
vations may or may not be a consequence of the glitch; 
the implied velocity of the disturbance, assuming that it 
originates near the pulsar, is ~ 0.7c, similar to the velocity 
inferred from changes in the Crab Nebula wisps. 

Future observations with Chandra and XMM can be 
used to gain further insight into the structure of the neu- 
tron star and its surrounding nebula. An observation 1.5 
yr after the glitch, now scheduled, will further constrain 
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models for the glitch and the parameters of the neutron 
star. Additional HRC observations will also be required 
to decide whether the nebula changes reported here are a 
consequence of the glitch, or whether they occur routinely 
in response to instabilities in the pulsar's relativistic wind, 
as appears to be the case in the Crab Nebula. Data from 
the EPIC PN camera on XMM will yield spectral clues 
helful in understanding the complex pulse profile, while 
either EPIC or Chandra's ACIS could be used to search 
for spectral changes caused by synchrotron energy losses 
and/or internal shocks in the nebula. 
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Fig. 1.- An X-ray image of the region containing the Vela pulsar PSR B0833 — 45. The two Chandra, High Resolution 
Camera (HRC) X-ray datasets (0.1-10 keV) centered on the pulsar were summed and smoothed with a Gaussian with 
a ~ 0.5'. The intensity scale is chosen to highlight the diffuse SNR emission surrounding the bright pulsar nebular which 
is fully saturated in this image. 
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The Vela Pulsar X-ray Nebula 



-45'^9'30"|= 



O 
O 

o 

+-> 

■1—1 

r— H 

Q 
Q 



-lO'O" - 



-10' 30 " - 



-U'O^' - 



-11 30"- 



■p-p-p 



■p-p-p 



■p-p-p 



■p-p 



■p-p-p 



■p-p 




J L 



J 



J \ L 



J l_L 



J \ L 



.^J L 



20" 



18' 



16^ 



Right Ascension (J2000) 



Fig. 2. — A close-up view of the region surrounding the Vela pulsar PSR B0833 — 45. The image includes the summed 
data from two epochs separated by a month and is centered on the pulsar and scaled to highlight the surrounding nebula 
emission. The data in this and the two subsequent figures have been smoothed with a Gaussian with a = 0.66". A 
toroidal structure and perpendicular jet similar to that seen in the Crab Nebula is apparent. Also evident is a faint halo 
of emission likely associated with the post-shock pulsar wind (see text). 
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Fig. 3. — A comparison of the relativistic wind nebulae surrounding two young pulsars observed by the Chandra 
Observatory: the 1000-yr Crab pulsar (right) and the 10 kyr Vela pulsar (left). The images are displayed with the 

same plate scale, although the Vela nebula is sixteen times smaller physically assuming distances of 2 kpc (Crab) and 250 
pc (Vela); the circles in the two images represent the same physical size at the respective pulsars. 
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Fig. 4. — The same image as in Fig. 2 with the various extraction regions described in the text shown: the circular 
extraction and background amilii region for the pulsar, the elliptical annulii bounding the toroidal structure which define 
the region used for the azimuthal profiles, the sector used to extract the radial profile (dashed lines), and finally, the 
background anulii (see text for details). The sector bounded by the the elliptical annulii and dashed lines is found to 
brighten between the two observations. 
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Fig. 5. — Top panel: the X-ray (0.1 — 10 keV) pulse profile of the Vela pulsar obtained with the Chandra HRC from 
the combined observations. Bottom panel: the difference profile - 2000 21 February minus 2000 20 January. There is no 
significant pulse shape change at any phase. The phase (^o (radio) is the location of the radio pulse. 
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Fig. 6. — Top panel: The radial profile of the Vela wind nebula in elliptical annuli. The profile is centered on the pulsar 

which has been suppressed and is restricted to the azimuthal sector delineated in Figure 4. Bottom panel: The difi^erence 
between the radial profiles for the two observations. Note the 7.8cr excess in the second observation between 13.5" and 
18" from the pulsar. 
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Fig. 7. — The azimuthal distribution of surface brightness in an elliptical annular ring between 13.5" and 18" from the 
pulsar (upper panel). The 7.6(t excess (see Fig 6) found in the second observation (lower panel) is located between position 
angles 240° and 10°. If caused by the glitch, this excess requires signal propogation at 0.7c. 




Fig. 8. — a) An idealized geometry for the Vela X-ray nebula as seen in the Chandra images. The two bright arcs are 
assumed to lie on the front surface of an equatorial, toroidal wind whose axis of symmetry is the rotation axis 17 of the 
neutron star. Lengths a and h are the semi-major and semi-minor axes of the projected ellipses, and s is the projected 
distance between them. For scaling, a = 25.7" translates to 0.96 x 10^^ cm at an assumed distance of 250 pc. The vector 
Vt is the direction of proper motion of the pulsar, b) A side view showing the definitions of the angles between the rotation 
axis fi, the magnetic axis /i, and the line of sight. The angle q — C is inferred from radio pulse polarization measurements, 
and ^ comes from deprojection of the X-ray structures. The angle a is only approximately known. 



